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We show that because of the multinucleon mechanism effects, the algorithm used to reconstruct 

^SJ ' the neutrino energy is not adequate when deaUng with quasielastic-like events, and a distortion of 

T-H , the total flux unfolded cross section shape is produced. This amounts to a redistribution of strength 

' from high to low energies, which gives rise to a sizable excess (deficit) of low (high) energy neutrinos. 

' This distortion of the shape leads to a good description of the MiniBooNE unfolded CCQE-like cross 

-H , sections published in Ref. However, these changes in the shape are artifacts of the unfolding 

r—{ • process that ignores multinucleon mechanisms. 

^ ; 
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^ ' I. INTRODUCTION 

Ph. 

O |. In most theoretical works the name Quasielastic (QE) scattering is used for processes where the gauge boson W is 
D ' absorbed by just one nucleon, which together with a lepton is emitted (see Fig. [Ija)). However, in the MiniBooNE 
7^ , . measurement of Ref. [l| , QE is related to processes in which only a muon is detected in the final state. Though this 
definition could make sense because ejected nucleons are not detected in that experiment, it includes multinucleon 
CNJ ' processes (see Fig. [Ijb))^ and others like pion production followed by absorption. However, it discards pions coming 
J> off the nucleus, since they will give rise to additional leptons after their decay (see Fig.[IJc)). The MiniBooNE analysis 
of the data corrects (through a Monte Carlo estimate) for some of these events, where in the neutrino interaction a real 
pion is produced, but it escapes detection because it is reabsorbed in the nucleus, leading to multinucleon emission. 

As firstly pointed out by M. Martini et al. [1, and corroborated by our group 0, the data of Ref. [l| 
correspond to the sum of the QE (absorption by just one nucleon), and the multinucleon contributions. For this 
■ reason, we will use the name QE-like to quote the MiniBooNE data of Ref. fT'|. Also for simplicity, we will often 
', refer to the multinucleon mechanism contributions, though they include effects beyond gauge boson absorption by 
I ' a nucleon pair, as 2p2h (two particle-hole) effects. The 2p2h contributions allows to describe Q the CCQE-like 
• • flux averaged double differential cross section da /dE^d cos 9^ measured by MiniBooNE with values of Ma (nucleon 
axial mass) around 1.03 ± 0.02 GeV that is usually quoted as the world average [l3|- This is re-assuring from the 
^> theoretical point of view and more satisfactory than the situation envisaged by some other works that described these 
}^ , CCQE-like data in terms of a larger value of Ma of around 1.3-1.4 GeV P, [lTI-[r3. 
_<^_- For the QE cross-sections, the predictions of the model that we employed in [5j agree quite well with those ob- 
tained/used in Refs. [1, 0, 0]) and both groups also agree on the relevant role played by the 2p2h mechanisms to 
describe the MiniBooNE data. We, however, differ considerably in the size (about a factor of two) of the multinucleon 
effects 01 • Thus, although Martini et al., predictions look consistent with MiniBooNE data, however our predictions, 
when the 2p2h contribution is included, would favor a global normalization scale of about 0.9 (see [B!]). This would 
be consistent with the MiniBooNE estimate of a total normalization error of 10.7%. In view of the disagreement, we 
should emphasized here that our evaluation in [1, Q , of these pionless multinucleon emission contributions to the cross 
section is fully microscopical and it contains terms, which were either not considered or only approximately taken into 
account in [3, 4, 6]. Indeed, the results of these latter works rely on some computation of the 2p2h mechanisms for 
the (e, e') inclusive reaction (iSj), whose results are used for neutrino induced processes. Thus, it is clear that these 
latter calculations do not contain any information on axial or axial-vector contributions. 



o 
in 



^ Note that the intermediate pion in this figure is virtual and it is part of the AA'' — )■ NN interaction inside of the nucleus. Indeed, one 
should consider a full interaction model for the in medium baryon-baryon interaction. Thus, for instance, the model of Ref. Q contains, 
besides pion exchange, p— exchange, and short and long range (RPA) correlations. 
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We would also like to point out that the simple phenomenological approach adopted in IS] to account for the 2p2h 
effects also reinforces the picture that emerges from the works of Refs. jSi, [6j]. Yet, a partial microscopical calculation 
of the 2p2h contributions to the CCQE cross section has been also presented in Refs. [l^ and [l^, for neutrino and 
antineutrino induced reactions, respectively. In these works, the contribution of the vector meson exchange currents in 
the 2p2h sector is added to the QE neutrino or antineutrino cross section predictions deduced from a phenomenological 
model (SuSA) [l7| based on the super-scaling behavior of electron scattering data. In [ll], and for the neutrino case, 
the SuSA-t-2p2h results were also compared with those obtained from a relativistic mean field approach. Although, all 
these schemes do not account for the axial part of the 2p2h effects yet, the preliminary results also corroborate that 
2p2h meson exchange currents play an important role in both CCQE neutrino and antineutrino scattering, and that 
they may help to resolve the controversy on the nucleon axial mass raised by the recent MiniBooNE data. This is not 
surprising, since these two-body currents, that arise from microscopic relativistic modeling performed for inclusive 
electron scattering reactions, are known to result in a significant increase in the vector-vector transverse response 
function in QE electron scattering data (T9l - [2l| . 

The study and comparison in detail of the different models used to describe 2p2h effects, though of great interest, 
is left for future research. We aim here at determining the possible influence of the 2p2h excitations on the process 
needed to extract neutrino energy unfolded cross sections from the measured flux-average data. 
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FIG. 1: Mechanisms for W absorption inside of a nucleus. 



Indeed, there still exists another feature of neutrino physics which deserves attention and that has motivated this 
work. Neutrino beams are not monochromatic. For QE-like events, only the charged lepton is observed and the only 
measurable quantities are then its direction^ (scattering angle 0^ with respect to the neutrino beam direction) and 
its energy i?^. The energy of the neutrino that has originated the event is unknown. Then, it is common to define a 
reconstructed neutrino energy E'lcc as. 



ME^ ~ ml/2 
M - E^, + \p^\cosB^ 



(1) 



which will correspond to the energy of a neutrino that emits a muon, of energy E^^ and three- momentum p^, and a 
gauge boson W that is being absorbed by a nucleon of mass M at rest. Namely, the usual reconstruction procedure 
assumes that we are dealing with a genuine quasielastic event on a nucleon at rest, ie. i?rec is determined by the 
QE-peak condition — — g^/2A/, where is the W four momentum. Note that each event contributing to the 
fiux averaged double differential cross section da/dE^dcos9^ defines unambiguously a value of i?icc- The actual 
("true") energy, E, of the neutrino that has produced the event will not be exactly i?rcc- Actually, for each E'roc, 
there exists a distribution of true neutrino energies that could give rise to events whose muon kinematics would lead 
to the given value of -Eroc. Several effects can infiuence this distribution. Firstly, the Fermi motion which broadens 
the QE peak and the Pauli blocking which cuts the low momentum response. These effects are well known'^, usually 
are under control and lead to very minor changes in the process of expressing observables as a function of the true 
neutrino energy. This is because genuine QE events produce true energy distributions quite narrow and strongly 



^ From now on, we will always identify the charged lepton with a muon. 

^ It is also common to consider some corrections in the definition of Erec in Eq. JTJ to account for the binding energy of the target nucleon 
in the nucleus, but these corrections turn out to be irrelevant for our discussion. 
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peaked around the expected £"100 values [22|. However multinucleon mechanisms, relevant for QE-like processes, can 
indeed distort the expected (QE-based) {E^ec,E) distributions, since they produce distributions quite flat and that 
do not peak around E'rec |22|. The effects of the inclusion of multinucleon processes on the energy reconstruction have 
been investigated in Ref. |22|. within their 2p2h model and also estimated in Ref. fsO^, using some simplified model 
for the multinucleon mechanisms. 

We will show in this work that 2p2h effects sizably distort the shape of the total CCQE-likc flux unfolded cross 
section, as a function of the neutrino energy. Indeed, we will see that these multinucleon mechanisms produce a 
redistribution of strength from high energy to low energies, which gives rise to a sizable enhancement of the number 
of events attributed to low energy neutrinos leading to a good description of the unfolded cross section given in 
However, we will show that these changes in the shape are artifacts of the unfolding process that ignores multinucleon 
mechanisms. 



II. EXCESS OF LOW ENERGY NEUTRINOS IN THE MINIBOONE CCQE-LIKE FLUX UNFOLDED 

CROSS SECTION DATA 

The QE+multinucleon mechanism model of Ref. @, with Af^ — 1.049 GeV, provides an excellent description of 
the MiniBooNE neutrino flux folded CCQE-like da /dT^d cos 6 f_i differential cross sections given in even when 
no parameters have been fitted to data, beyond a global scale, A ^ 0.9. This scale A is consistent with the global 
normalization uncertainty of around 10% acknowledged iii ll[. 

The QE contribution used in Q was derived in Ref. (23l |T and it incorporates several nuclear effects. The main 
one is the medium polarization (RPA), including A-hole degrees of freedom and explicit tt and p meson exchanges 
in the vector-isovector channel of the effective nucleon-nucleon interaction. The model for multinucleon mechanisms 
has been fully discussed in Ref. ^"^^ it is based on a model for neutrino pion production derived in Refs. [2^. [25|. 
The whole model constitutes a natural extension of previous studies of photon, electron, and pion interactions with 
nuclei 

The prediction of the model, used in Ref. [sj, for the total flux unfolded neutrino CCQE-like cross section is depicted 
in Fig. [21 Several remarks are in order here: 

• The 2p2h contributions clearly improve the description of the data in Fig. [21 which are totally missed by the 
QE prediction. Though the model provides a reasonable description, we observe a sizable excess of low energy 
neutrinos in the data, that is not even covered by the theoretical error band. 

• As discussed above, the flux folded double differential cross-section data is well described in Ref. @, except for 
the global scale, A 0.9, which needs to be introduced there. It is to say, predicted cross sections in Ref. Q are 
globally around 10% smaller than the measured ones. This disagreement could be due to a theoretical under- 
estimation of the absolute number of neutrinos in the MiniBooNE flux. Thus, we should expect our predictions 
for the total unfolded cross section to under-cstimate the data points by about 10%, as well. However, we do 
not see this in Fig. [21 There is a problem in the neutrino-energy shape, as pointed out above, which would not 
be improved by increasing the size our predictions by a global factor. 

In any case, given that the actually measured quantity it is the double differential cross-section and that 
observable is well described by our model, any difference on the unfolded cross section must come from the 
unfolding procedure. 

• Finally, we should mention that the QE theoretical results for the cross section shown here (model from [51]) 
slightly differ from those in Fig. 18 (left) of our previous work of Ref. 0. The main difference is the inclusion of 
relativistic corrections. We discuss this point in some detail in Appendix[3 since we do not want to deviate here 
the attention from the main point of this work: because of the 2p2h effects, the algorithm used to reconstruct 
the neutrino energy is not adequate when dealing with QE-like events, and that it produces a distortion of the 
total CCQE-like flux unfolded cross section shape. 

Nevertheless, we should mention here that in Ref. [H and to account for Final State Interactions (FSI), we used 
the non-relativistic QE model of Ref. [2^, while in Fig. [21 the results depicted are those obtained from the 
relativistic model of [23] for the QE process, without the inclusion of FSI effects'*. This improvement in the 
model to account for the relativistic effects is the reason for the differences mentioned above. As a flnal remark. 



4 In Ref. 0, FSI effects are being treated within the non-relativistic scheme derived in Ref. [30(1 ■ A non-relativistic treatment is unsuitable 
for the large momenta transferred that are reached in the MiniBooNE neutrino flux folded da /dT^d cos 8 differential cross sections, but 
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FIG. 2: Flux-unfolded MiniBooNE v^^ CCQE-like cross section per neutron as a function of neutrino energy (data points) 
from Ref. Jj], together with the predictions derived from the model used in Ref. 5]. The yellow band accounts for theoretical 
uncertainties, as discussed in Ref. 9|, while the QE contribution includes relativistic effects, and some nuclear corrections, 
among other those due to long range RPA correlations. 



we stress that the results of Ref. Q for the double differential cross section and those displayed in Fig. [3] have 
been calculated with the same model. 

We see in Fig. [2] that the proportion of multinucleon events contributing to the QE-like signal is quite large in the 
whole energy range relevant in the MiniBooNE experiment. This questions the validity of the algorithm used to 
reconstruct the neutrino energy in Eq. ([l} . In the next section, we will explore in detail this problem and we will find 
out the shape distortion effects induced by the use of Eq. ([J) in the MiniBooNE data. 



III. THE DISTRIBUTION OF RECONSTRUCTED ENERGIES IN CCQE-LIKE EXPERIMENTS 

Let P(E!^^^E) be the conditional probability density of measuring an event with reconstructed energy (combination 
of muon energy and scattering angle given in Eq. (jlj) comprised in the interval [£'rcci -^rcc + "^-^rcc] and induced by the 
interaction with the nuclear target of a neutrino of energy E (it is to say, conditional probability density of obtaining 
i?rec "given" E). This probability density can be computed theoretically as: 

P{EUE) = J_^(£;;iJ_ = ijo^j (2) 
a{E) dErcc 

where (t{E) is the integrated CCQE-like cross section for neutrinos of energy E, and the distribution da/dE^-cc is 
obtained from the double differential cross section, with respect the energy and scattering angle of the outgoing 
muon, as 



da _ „n . /"^ ,„ rfV 



dErcc Jra„ dE^ccdE 



(i?;£;L)= dE,———{E;E^J^ dE, 



9(cos e^,) 



dE„ 



d(cos 9^j)dE^ 



{E;E^J (3) 



for a fix value of the reconstructed energy E^cc = E^^^ and "true" neutrino energy E. Eq. ([T|) can be used to express 
cos in terms of -E^i and -Ercc- Besides, the Jacobian can be trivially computed also from Eq. ^ and it reads 

d{cos0f,) ^ ME^ - ml/2 
dErcc EUp^\ ^ > 



it is more appropriated for the total unfolded cross section as long as the neutrino energy is sufficiently small. A final consideration, in 
Ref. [H, it was found that the main effect of FSI is a shift of ~ 10 MeV of the QE peak for neutrino energies closer to the MiniBooNE 
neutrino flux mean energy, (E) ~ 800 MeV and that it has little impact on the integrated cross section (we will illustrate, within the 
model of Ref. [23| . this latter affirmation also in the Appendix |Aj . Ifowever, we cannot discard the possibility that these effects could 
be more important in the angle and energy distributions for low energy neutrinos. Moreover, it has been also pointed out that some 
relativistic approaches to account for FSI lead to larger variations of the total cross section |31| . 
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We would like to stress that da / dEicdE; E^^^) is an observable, but this distribution is not accessible to experiments 
where only the kinematics of the outgoing muon is measured. 

On the other hand, let Pi.ec(£'roc) be the probability density of measuring an event with reconstructed energy E^cc, 

Proc(Srcc) = J P{E,,,\E)Ptr-nc{E)dE (5) 

where 

Ptruc(£;) = -j^<PiE)a{E), (a) = / 'i>iE')aiE')dE' (6) 

is the density probability of having an event due to the interaction of a neutrino, with energy between E and E + dE, 
with the nuclear target. $ is the neutrino flux normalized to one, and (cr) is the total flux averaged cross section. It 
trivially follows, 

1 f dcT 

Proc{E,,c) = rr / :77^(^^;^rec)*(£;)d-B (7) 

[a) J ai^rcc 

a magnitude that can be measured in a CCQE-like experiment. 

IV. THE DISTRIBUTION OF TRUE NEUTRINO ENERGIES Ptruc(£;) FROM CCQE-LIKE 

EXPERIMENTS 

We use Bayes's theorem to estimate Ptruo(-£') from the measured density probability Prcc(-£'rcc)- To that end, let us 
introduce P{E\E,-qc) that is, given an event of reconstructed energy -Eiccj the conditional density probability of being 
produced by a neutrino of energy E. It follows, 

P„^,iE>^l,EME,,.)P,ElE.,.> (8) 

Now, since Bayes's theorem reads 

PiE\E,.^^) = ^(^.^c|i^)Pt.Uc(i?) 

we deduce 

P(F\p , ^E)da/dE,,,{E;E,,,) 

n^\J^roc) - jdE"<f{E")da/dEUE";E,,,) ^'"^ 

from Eqs. ([U, ([S]) and ([7]). The recent work of M. Martini et al. [22| pays an special attention to this distribution 
that, as we observe, depends on the neutrino flux. The above equation implies 

P rp^- f^P P (P ^ HE)da/dE,,,{E;E,,,) 

nruo[J^) - J CXAec/^rccl^-rcc j j ^E" ^{E")da / dE,,,{E" ■ E,,,) ^ ' 

Finally and attending to the existing relation between Ptiuc{E) and (j{E) in Eq. (jG)), we could write 

da/dE,,,{E;E 

roc J 



a{E) = / dE,,, {a)P,,,{E,,,) 



(12) 



/ dE"^E")da/dE,,,{E"-E,,,) 
A consistency check is obtained if we substitute Eq. (O in Eq. (fT2)) which leads to 

dE,,,—^{E;E,,c) (13) 

that it is trivially satisfied thanks to the definition of da / dE^cc{E\ i?rec) in Eq. ([3]). 

Eq. (fT2l) might be used to estimate the integrated flux unfolded cross section from data. CCQE-like experiments 
measure the quantities that appear in the first bracket of this equation, namely, (a) x Prec(^rec)- We have already 
discussed about Pi.oc(-Eroc), while the total flux averaged cross section (cr) is determined by the ratio of the total 
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number of events {Nevent) over the number of incident neutrinos per unit of area {Nine)- Nevent is directly measured 
while for Ni^c there exist, in principle, accurate theoretical predictions. Note that the flux $(£'), normalized to one, 
gives only the shape of the neutrino flux, but it is independent of the total number of incident neutrinos. 

Thus, if one had a theoretical model for the second bracket [dcr/(ii?rcc(£'; iSicc)/ / dE"^{E")da/dEi-cc{E";E^cc)] 
in Eq. ([T2|) or equivalently for P{E\Ercc), one could extract the flux unfolded cross section (t{E) after folding it 
with the measured data ((tr) x Pi.cc(i?icc))- This method reproduces [sll the data unfolding used by the MiniBooNE 
collaboration in Ref. H and described in Ref. f36l|. The iterative unfolding method in MiniBooNE is needed due to 
the statistical fluctuations in the data and the lack of "a priori" knowledge of the Ptrue{E) probability in Eq. ([S]), 
both are not relevant for theoretical calculations. However, as a proof of the validity of this approach we have checked 
that the iterative method in [36i| yields to identical results. As discussed in the introduction, the works of Refs. 0-11] 
show that the QE-like (and/or differential) cross section is given by the sum 



cj{E)^a'^^{E) + a''^'''\E) 



2p2h/ 



(14) 



of the genuine QE and the multinucleon contributions. Up to now, experimental analysis have completely neglected the 
latter (2p2h) cross section, while well established nuclear corrections, like RPA correlations, have also been ignored 
in the computation of the former one (QE). As a consequence, a high value of Ma > 1.3 GeV is required in the 
MiniBooNE analysis to describe the flux folded da/dq^ or da /dT^dcosO^ distributions [1]. But, if this approximate 
model is used in Eq. (|T2|) to extract the unfolded cross section. 



o-appx(i^) = / dE^cc {a)Prcc{E, 



) X 

J Exp 



da/dE^cc{E; E^, 



J dE"<i>{E")da/dE,UE"\E,,,) 



(15) 



QE no RPA, A/a>1-3 GcV 



the resulting estimate (Tappx(-E') might signiflcantly differ from the real QE-like neutrino-nucleus cross section cj{E). 
Actually, we will show that a redistribution of strength from high energy to low neutrino energies is being produced. 
To illustrate this, we will focus on the total cross section data on carbon published by the MiniBooNE in |l|. We take 
in Eq. ([T5|) 



(CT)^'roc(-E^rcc) 



Exp 



/ da 



V dE,c 



{E';E,, 



A/a = 1.049 GcV 
QE+RPA, 



dE,, 



{E';Eroc))HE')dE' 



(16) 



where we have used Eq. ([7]) with our best theoretical model, since we should mimic the experiment. Indeed, this 
is the same model we employed in Ref. Q to successfully describe the MiniBooNE CCQE-like flux averaged double 
differential cross section da /dE^dcosO^^ up to a global normalization scale A (= 0.89 ± 0.01). This latter parameter 
is the only one which is fltted to data (x^/do/ — 53/137). There are no free parameters in the description of nuclear 
effects, since they were flxed in previous studies of photon, electron, and pion interactions with nuclei [26l-l29i ISTllssj . 
Besides, form factors are determined in independent analysis of the experimental data on nucleons. In particular, 
the model uses a value for the nucleon axial mass of Ma — 1.049 GeV, that agrees within errors with the world 
average [1, [l^ value of 1.03 ± 0.02 GeV. As already mentioned, the genuine QE piece was computed in Ref. [2], it 
includes relativistic effects, and some other nuclear corrections in addition to the RPA ones, explicitly included in 
the label of Eq. p^ . though it does not include FSI effects. The multinucleon cross section is taken from 0. The 
QE and 2p2h contributions to the integrated cross section a{E) in carbon were displayed in Fig.[5J For simplicity, in 
what follows, wc will label this QE model as "QE (rel-fRPA)" as in that figure. 

On the other hand, to compute the second factor in the right hand side of Eq. (|15p we need to mimic the input 
used in the experimental analysis. To that end, we use a simple Fermi gas model with Ma = 1.32 GeV that only 
accounts for the genuine QE contribution and that does not include RPA corrections. The value of Ma in this model 
was fitted to the flux-folded double differential cross section da/dE^dcos6^ in Ref. Q, leading to an excellent value 
of /dof = 35/137. This model should be quite similar to the one originally used in the MiniBooNE analysis. The 
main difference being that we use a local rather than global Fermi gas in the calculation. Gathering the different 
terms, we have 



QE (rcl+RPA) 
appx 



{E) 



2p2h 
appx 



[E) 



(17) 



QE (rcl+RPA) 
•^appx 



{E) 



dEr, 



A/a = 1.049 GcV 



^^QE (rcl+RPA) 



{E';E,,,)^E')dE' X 



QE no RPA, Ma = 1.32 GcV 



da/dE^cc{E; E'rcc) 



/ dE"<i>{E")da/dE,,,{E"-E,,,) 



(18) 
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FIG. 3: Theoretical (a) and approximate (o"appx, defined in Eqs. (I17|l - (|19p ) CCQE-like integrated cross sections in carbon as 
a function of the neutrino energy. 





FIG. 4: Differential cross section da /dErcc{E; E^cc) as a function of the "true" neutrino energy (E) and four different values of 
Eiec = 0.5, 0.75, 1 and 1.25 GeV, which are indicated by the vertical lines. QE contributions, from the (rel+RPA) model, are 
displayed in top panel, while 2p2h ones are shown in the bottom plot. 



crl^l^iE) - J dE,,,J ^—{E';E,,,ME')dE' x 

QE no RPA, Ma = 1-32 GcV 

Results both for the QE and the 2p2h contributions to CTappx are shown in Fig. [31 As can be appreciated there, 
o'uppxiE) is an excellent approximation to the real (7(E) cross section in the case of the QE contribution. The reason 
is that for genuine QE processes, the distribution da'^^ / dEicc{E; Eycc) is strongly peaked around E w E'roc, as seen 



da/dE,,,{E;E,, 



r dE"<i>{E")d(T/dE,^^{E"; E„ 



(19) 
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FIG. 5: Theoretical a and approximate (Tappx, ( defined in Eqs. (|17p - (|19|l l CCQE-like integrated cross sections in carbon as a 
function of the neutrino energy. For consistency with our previous results for the flux-folded double differential cross section 
da / dE^dcos 6^ in Ref. the MiniBooNE data [ij and errors have been re-scaled by a factor 0.89. The shape errors are shown 
for the MiniBooNE data. 



in the top panel^ of Fig. Sl^ 

Though the actual da^^ / dE^cc distributions have some widths, these differential cross sections are sufficiently 

peaked to render the width effects on the ratio (^a'^^x'^^^^^^^^ / (rGi-i-RPA)~j quantitatively irrelevant, as the results 

of the Fig. [3] indicate. Thus, we could conclude that when dealing only with genuine QE events the procedure outlined 
in Eq. to obtain the flux unfolded cross section is quite accurate. This is despite of the fact that RPA correlations 
and other nuclear effects were not considered in the ansatz for P{E\Ercc) (second bracket in Eq. (|12l) '). Note however, 

that all nuclear effects are included in the first factor (a) Pi-cc{Ei-cc) in Eq. 

L J Exp 

However, the situation is drastically different for the 2p2h contribution case, as one can also observe in Fig. |3l 
Indeed, it turns out that a'^^p^{E) is a poor estimate of the actual multinucleon mechanism contribution (7^P^^(i?). 
As before, if we approximate 



^ Actually the peak is shifted about 25 MeV up to higher energies since the bound energy of the target nucleon is not considered in Eq. JT}. 

On the other hand, RPA correlations modify the size but do not affect significantly this peak structure. 
^ If for illustration purposes, we use a Dirac's delta to approximate ^ — {E; E^cc) ~ <t*^^(_E) x 5{E — E^cc) then, we will have 



QE no RPA, Ma = 1.32 GcV 



da/dErcciE; E^cc) 



JdE"'!>{E")da/dErcc{E"-E,cc) 
and therefore, independently of the nuclear model for QE, within this limit 

A/_4 = 1.049 GcV 

5{E - Ere, 



S{E - Erec) 
^Eroc) 



QE (rel+RPA) , , 
-'^appx \^ } ■ 



dEr, 



^(Ercc) 



^QE (rcl + RPA) 



(20) 



(21) 
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Taking into account that da'^^'^^ /dEycc{E; -Ercc) is almost negligible for E < E^^c and that it shows a quite long tail 
above this energy (see bottom panel of Fig. 0]), is then easy to understand the redistribution of strength from high 
to low neutrino energies observed in a'^^{E) when it is compared to the actual (T^P^'^(i^) cross section. Up to some 
approximation the area is conserved, though. 

Finally, in Fig.[5]we compare the MiniBooNE CCQE-like data with both a and Cappx. We see an excellent agreement 
between the latter one and the data scaled down by a factor 0.89. As mentioned, our QE(rcl+RPA)+2p2h model 
successfully describes the MiniBooNE CCQE-like flux averaged double differential cross section da /dE^d cos 9^ data, 
up to a global scale A (= 0.89 ± 0.01) This is the best observable to compare with theoretical models because 
both the muon angle and energy are directly measured quantities, and thus the shape of this distribution is readily 
obtained from the number of events measured for each muon kinematical bin. To obtain the absolute normalization 
of the distribution, however, it is necessary to rely on some estimate for the number of incident neutrinos per unit of 
area (Nine)- We believe the obtained value for A ~ 0.89 in Q indicates that the actual number of incident neutrinos 
per unit of area might be larger than the central value assumed in the MiniBooNE analysis. This would be still 
consistent with the MiniBooNE estimate of a total normalization error of 10.7% [H. The value of Nine is needed to 
estimate (cr) Exp in the expression of Eq. (|15|) for crappx(£^), and thus the predictions of our model in Eq. ()16|) would 
have to be multiplied by 1/A, or equivalently the data have to be scaled down by a factor A. 

Coming back to the results displayed in Fig. [5l we should conclude that the unfolded cross section published in [l[ 
appreciably differs from the real one cf{E). Actually, it is not a very clean observable after noticing the importance of 
multinucleon mechanisms, because the unfolding itself is model dependent and assumes that the events are purely QE. 
The same limitation occurs for the differential cross section da/dq^, given that is also deduced assuming the events 
are QE. When compared with the "real" <t{E), the MiniBooNE unfolded cross section exhibits an excess (deficit) of 
low (high) energy neutrinos, which is an artifact of the unfolding process that ignores multinucleon mechanisms. 

The semi-phenomenological model of Refs. [1, 0] predicts a theoretical cross section a{E) that provides a good 
description of the MiniBooNE unfolded cross section of Ref. [l|. However, we have shown here that these data do 
not correspond to the actual cross section because the unfolding process is biased. To compare with these data, the 
authors of [3, ^ should carried out a procedure similar to that proposed in Eq. (|15p . Since the model of these works 
includes strong 2p2h contributions, we would expect an appreciable change in the shape, as discussed above, that 
might jlistort quantitatively and qualitatively the agreement with the MiniBooNE unfolded cross section found in 
Refs. a 141. 

In Ref. [22, Martini et al. have paid special attention to the flux dependent P{E\E,-cc) probability. We however 
believe that P{E^ec\E) (or equivalently the differential cross section da/dE^cc) could be more illuminating. First, 
because it does not depend on the neutrino flux, and second because it can be used by experiments to determine <t{E) 
thanks to Bayes's theorem. Finally, we should also mention a recent and quite comprehensive work on neutrino-nucleus 
observables [39| has also found, using some simple models for multinucleon mechanisms, that 2p2h interactions lead 
to a downward shift of the reconstructed energy in agreement with our results. 



V. CONCLUSIONS 



We have shown that because of the the multinucleon mechanism effects, the algorithm used to reconstruct the 
neutrino energy is not adequate when dealing with quasielastic-like events. This effect is relevant to neutrino oscil- 
lation experiments that uses the CCQE samples to compute the neutrino energy. The assumption of a pure CCQE 
interaction introduces biases in the determination of Am^ and mixing angle. Moreover, 2p2h contributions put also 
limitations on the validity of the flux unfolding procedure used in . The MiniBooNE unfolded cross section exhibits 
an excess (deficit) of low (high) energy neutrinos, which is an artifact of the unfolding process that ignores multin- 
ucleon mechanisms. Actually, Cappx, defined in Eqs. PT)) - (ITO1) . provides an excellent description of the data of [l|. 
This, together with our previous results in Ref. for the CCQE-like fiux averaged double differential cross section 
da / dE ^d cos 9 fj,, make us quite confident on the reliability of our QE(rel-|-RPA)+2p2h microscopical model derived 
in Refs. |2, H^. Furthermore, because it is just a natural extension of previous successful studies of photon, electron, 
and pion interactions with nuclei [26l - [29j . 



Appendix A: Relativistic vs non-relativistic QE total CC neutrino cross sections within the scheme of 

Ref. [23| 

The left panel of Fig. 18 in Ref. 0] corresponds to the flux-unfolded MiniBooNE CCQE-like cross section per 
neutron as a function of the neutrino energy. There, both the QE and the 2p2h contributions to the total cross section 
were also shown separately. The latter ones were computed fully relativistically, while the QE predictions were taken 
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FIG. 6: Different theoretical predictions for neutrino CCQE total cross section off ^^C obtained from the model of Ref. [2^ . 
Yellow bands account for a 15% theoretical uncertainties that might affect the nuclear corrections included in the model of 
Ref. [l^], as discussed in [ssl ]. 



from a previous work [23|. Concretely, results that included RPA and FSI effects, were selected and displayed in 
the Fig. 18. Since the approach used in [2^ to account for FSI effects was not relativistic, the QE curves displayed 
in (both panels) Fig. 18 of Ref. Q neglect some relativistic effects for the nucleons. In particular, those results 
are based on a non-relativistic approximation for the nucleon (particle and hole) propagators''. This is one of the 
sources of systematic errors, among others, that should be accounted for by the bands of theoretical uncertainties 
displayed in Fig. 18 (see the discussion of the third paragraph of pag. 16 in [2]). The use of non-relativistic nucleon 
propagators is responsible of the odd behaviour of the QE results in the Fig. 18 of Ref. when the neutrino energy 
increases. Indeed, the QE cross sections are too big for neutrino energies above 0.5-0.6 GeV and these predictions 
clearly depart 33] from the common pattern exhibited by different models collected in a review talk presented in the 
NUINT 2009 Workshop [H. 

In the top panel of Fig. [6l we show the size of FSI and relativistic effects within the QE model of Ref. 23] (note 
that there, results obtained using relativistic nucleon propagators, but neglecting FSI, were also shown). In all curves 
of this top panel, RPA effects are taken into account. We see that FSI have little effect on the integrated cross 
sections, though FSI might affect differential distributions [l^, and we should also pointed out that some relativistic 
approaches to FSI lead to larger increases of the total cross section ^sTj. On the other hand, for neutrino energies 
above 1 GeV, relativistic effects for the nucleons reduce the cross sections by around 15%, but still the relativistic QE 



Actually, all QE results of Ref. were obtained with non-relativistic nucleon propagators. 
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results lie within the theoretical error band assumed for the QE predictions in Ref. 

The flux folded CC double differential neutrino cross section, measured by the MiniBooNE Collaboration, was 
analyzed in Ref. Q by using the full relativistic model of Ref. [2^ without the inclusion of FSI, but taking into 
account RPA correlations, which effects in the integrated flux-unfolded cross section can be seen in the bottom panel 
of Fig. ini This is the model that is used in this work. Note, as can be appreciated in this latter plot, that though 
RPA effects are quite relevant for low neutrino energies, they are negligible above 1 GeV, and much smaller than the 
theoretical uncertainties, discussed in [s^, for E > 0.7 GeV. 



Acknowledgments 



We thank R. Tayloe and G. Zeller for useful comments regarding the data unfolding used by the MiniBooNE 
coUaboration. This research was supported by DGI and FEDER funds, under contracts FIS2011-28853-C02-01, 
F1S2011-28853-C02-02, FIS2011-24149, FPA2011-29823-C02-02 and the Spanish Consolider-Ingenio 2010 Programme 
CPAN (CSD2007-00042), by Generalitat Valenciana under contract PROMETEO/2009/0090, by Junta de Andalucfa 
grant FQM-225 and by the EU HadronPhysics2 project, grant agreement no. 227431. 



A. A. Aguilar-Arevalo et al. [MiniBooNE Collaboration], Phys. Rev. D 81 (2010) 092005. 

J. Nieves, I. Ruiz Simo and M. J. Vicente Vacas, Phys. Rev. C 83, 045501 (2011). 

M. Martini, M. Ericson, G. Chanfray, J. Marteau, Phys. Rev. C80 (2009) 065501. 

M. Martini, M. Ericson, G. Chanfray, J. Marteau, Phys. Rev. C81 (2010) 045502. 

J. Nieves, I. R. Simo and M. J. V. Vacas, Phys. Lett. B 707, 72 (2012). 

M. Martini, M. Ericson and G. Chanfray, Phys. Rev. C 84, 055502 (2011). 

J. Nieves talk at the Workshop 'Modeling Charge-changing and Neutral-current Neutrino Reactions with Nuclei', ECT*, 
Trento, Italy, december 2011. 

W. M. Alberico, M. Ericson and A. Molinari, Annals Phys. 154, 356 (1984). 
V. Bernard, L. Elouadrhiri and U. .G. Meissner, J. Phys. G G 28, Rl (2002). 
V. Lyubushkin et al. [NOMAD Collaboration], Eur. Phys. J. C 63, 355 (2009). 
O. Benhar, P. Coletti and D. Meloni, Phys. Rev. Lett. 105 (2010) 132301. 

C. Juszczak, J. T. Sobczyk and J. Zmuda, Phys. Rev. C 82 (2010) 045502. 
A. V. Butkevich, Phys. Rev. C 82 (2010) 055501. 

O. Lalakulich, K. Gallmeister and U. Mosel, arXiv: 1203.2935 [nucl-th]. 

J. E. Amaro, M. B. Barbaro, J. A. Caballero, T. W. Donnelly and C. F. Williamson, Phys. Lett. B 696 (2011) 151. 

J. E. Amaro, M. B. Barbaro, J. A. Caballero and T. W. Donnelly, Phys. Rev. Lett. 108 (2012) 152501. 

J. E. Amaro, M. B. Barbaro, J. A. Caballero, T. W. Donnelly, A. Molinari and I. Sick, Phys. Rev. C 71 (2005) 015501. 

J. E. Amaro, M. B. Barbaro, J. A. Caballero, T. W. Donnelly and J. M. Udias, Phys. Rev. D 84 (2011) 033004. 

J. E. Amaro, C. Maieron, M. B. Barbaro, J. A. Caballero and T. W. Donnelly, Phys. Rev. C 82 (2010) 044601. 

J. E. Amaro, M. B. Barbaro, J. A. Caballero, T. W. Donnelly, C. Maieron and J. M. Udias, Phys. Rev. C 81 (2010) 

014606. 

A. De Pace, M. Nardi, W. M. Alberico, T. W. Donnelly and A. Mohnari, Nucl. Phys. A 726 (2003) 303. 
M. Martini, M. Ericson and G. Chanfray, Phys. Rev. D 85, 093012 (2012). 

J. Nieves, J. E. Amaro and M. Valverde, Phys. Rev. C 70 (2004) 055503 [Erratum-ibid. C 72 (2005) 019902]. 
E. Hernandez, J. Nieves and M. Valverde, Phys. Rev. D 76 (2007) 033005. 

E. Hernandez, J. Nieves, M. Valverde and M. J. Vicente Vacas, Phys. Rev. D 81 (2010) 085046. 

A. Gil, J. Nieves and E. Oset, Nucl. Phys. A 627 (1997) 543. 

R. C. Carrasco and E. Oset, Nucl. Phys. A 536 (1992) 445. 

J. Nieves, E. Oset and C. Garcia-Recio, Nucl. Phys. A 554 (1993) 509. 

J. Nieves, E. Oset and C. Garcia-Recio, Nucl. Phys. A 554 (1993) 554. 

P. Fernandez de Cordoba and E. Oset, Phys. Rev. C 46 (1992) 1697. 

A. Meucci, M. B. Barbaro, J. A. Caballero, C. Giusti and J. M. Udias, Phys. Rev. Lett. 107 (2011) 172501. 

R. Tayloe, T. Katori and G. Zeller, private communication. 

We are very much indebted to J. A. Caballero who pointed out this fact. 

S. Boyd, S. Dytman, E. Hernandez, J. Sobczyk and R. Tacik, AIP Conf. Proc. 1189 (2009) 60. 

M. Valverde, J. E. Amaro and J. Nieves, Phys. Lett. B 638 (2006) 325. 

G.D'Agostini, Nucl. Inst. Meth. A 362(1995) 487. 

E. Oset, H. Toki and W. Weise, Phys. Rept. 83 (1982) 281. 

L. L. Salcedo, E. Oset, M. J. Vicente- Vacas and C. Garcia-Recio, Nucl. Phys. A 484 (1988) 557. 
U. Mosel and O. Lalakulich, arXiv: 1204.2269 [nucl-th]. 

D. Meloni and M. Martini, arXiv: 1203.3335 [hep-ph]. 



